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ABSTRACT 

We present the results of Chandra and Suzaku X-ray observations of nine moderate-redshift 
(0.16 < z < 0.42) clusters discovered via the Red-sequence Cluster Survey (RCS). Surface 
brightness profiles are fitted to beta models, gas masses are determined, integrated spectra 
are extracted within R2500, and X-ray temperatures and luminosities are inferred. The L x - 
Tx relationship expected from self-similar evolution is tested by comparing this sample to our 
previous X-ray investigation of nine high-redshift (0.6 < z < 1.0) optically selected clusters. 
We find that optically selected clusters are systematically less luminous than X-ray selected 
clusters of similar X-ray temperature at both moderate and high-z. We are unable to constrain 
evolution in the Lx-Tx relation with these data, but find it consistent with no evolution, within 
relatively large uncertainties. To investigate selection effects, we compare the X-ray properties 
of our sample to those of clusters in the representative X-ray selected REXCESS sample, also 
determined within R2soo- We find that while RCS cluster X-ray properties span the entire 
range of those of massive clusters selected by other methods, their average X-ray properties 
are most similar to those of dynamically disturbed X-ray selected clusters. This similarity 
suggests that the true cluster distribution might contain a higher fraction of disturbed objects 
than are typically detected in X-ray selected surveys. 

Key words: cosmology: observations — cosmology: large scale structure — X-rays: galax- 
ies: clusters — galaxies: clusters: general. 



1 INTRODUCTION 

By virtue of their size, clusters of galaxies are an important 
source of information about the underlying cosmology of the uni- 
verse (e.g., Allen et al. 2011, and references therein). Since the 
advent of the first cluster survey over 50 years ago (Abell 1958), 
the search for these impressive objects has grown in cosmologi- 
cal impact, and is currently being pursued with the ultimate ob- 
jective of constraining w, the dark energy equation of state. Im- 
proved technology has enabled large-area cluster searches in sev- 
eral wavebands (e.g., X-ray, optical, millimetre, submm, radio) and 
out to formative redshifts. Consequently, the number of recent, un- 
derway, and planned cluster surveys is staggering: Bohringer et 
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al. 2004 (REFLEX); Valtchanov et al. 2004 (XMM-LSS); Glad- 
ders & Yee 2005 (RCS); Wester & Dark Energy Survey Collab- 
oration 2005 (DES); Content et al. 2008 (EUCLID); Wilson et 
al. 2009 (SpARCS); Conconi et al. 2010 (WFXT); Vanderlinde et 
al. 2010 (SPT); Lloyd-Davies et al. 2011 (XCS); Marriage et al. 
2011 (ACT); Planck Collaboration VIII 2011 (Planck); Predehl et 
al. 2011 (eROSITA); Suhada et al. 2012 (XMM-BCS) and many 
more. 

One overarching goal of these endeavors is to chart the evolu- 
tion of the cluster mass function, thereby providing key constraints 
on the progression of large-scale structure formation in the uni- 
verse. Constructing a mass function requires two key elements: the 
ability to find clusters, and an efficient method of mass estima- 
tion. Obtaining either piece necessitates an accurate understanding 



© 0000 RAS 



2 A. K. Hicks et al. 



Table 1. Cluster Sample 
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a Spectroscopic (Blindert 2006) 

^ From X-ray spectra (this work), see text. 

c Spectroscopic (Ellingson et al. 2012) 

d Photometric 



of relationships between the observable properties of clusters (i.e., 
baryons) and their underlying dark matter distributions. 

The baryonic mass in clusters is dominated by the hot intra- 
cluster medium (ICM), therefore some of the most powerful meth- 
ods of cluster identification (X-ray, SZ) use the ICM to both find 
clusters and estimate their masses. Several processes occur within 
clusters, however, that can alter the physical characteristics and 
distribution of cluster baryons. These include feedback (e.g., Voit 
& Donahue 2005), mergers (e.g., Chatzikos 2012), and radiative 
cooling (e.g., Pratt et al. 2009), among other non-gravitational pro- 
cesses (e.g., Nagai 2006). Additionally, the frequency with which 
those processes occur may vary with redshift (e.g., Barger et al. 
2005). 

Understanding the gamut of ICM properties requires X-ray 
observations of samples chosen independently of their X-ray char- 
acteristics. This has indeed been pursued multiple times (e.g., 
Holden et al. 1997; Donahue et al. 2001, 2002; Basilakos et al. 
2004; Gilbank et al. 2004; Popesso et al. 2004; Sanchez et al. 2005). 
Most of these surveys have identified a significant population of 
low-Lx clusters. Unfortunately, there is a striking dearth of ade- 
quate X-ray data for these samples. Without convincing L x and 
T x measurements it is impossible to verify alternatively obtained 
cluster masses or investigate trends in core gas density. Due largely 
to this lack of quality X-ray data, concrete physical explanations 
for the observed scatter in cluster properties have not yet been well 
determined. 

In our previous work (Hicks et al. 2008) we use pointed 
Chandra observations to compare the ICM properties of high-z op- 
tically selected clusters (RCS; Gladders & Yee 2005) with those 
of moderate-redshift X-ray selected clusters (CNOC; e.g., Yee, 
Ellingson, & Catiberg 1996). Our results indicate that the typical 
central (~R25ool where p c imt/Pcnt = 2500) gas mass fractions of the 
RCS sample are strikingly lower than those found in X-ray selected 
systems. This is also illustrated by discrepancies in the normal- 
ization of the L x -T x relationship between the two samples. These 
comparisons, however, were made between clusters of differing av- 
erage redshift (zcnoc ~ 0.3, zrcs ~ 0.8), which were gathered via 
significantly different means of selection, making it difficult to dis- 
entangle the effects of selection bias from possible signatures of 
cluster evolution. 

The main objectives of this work are to isolate these poten- 
tial factors from one another and identify the cause of discrepan- 



cies in ICM properties between optically and X-ray selected cluster 
samples. In Section 2 we introduce the present cluster sample and 
describe the initial processing of our Chandra and Suzaku X-ray 
observations. In Sections 3 and 4 we probe the surface brightness, 
temperature, metallicity, and density of the hot ICM present in each 
cluster. In Section 5 we compare the ICM properties of RCS clus- 
ters to X-ray selected samples, including the R EXCESS sample 
for which quantities have been extracted within R25oo- Finally, we 
compare X-ray temperature to velocity dispersion for a subset of 
our targets (Section 6). Our results are summarized in Section 7. 

Unless otherwise noted, this paper assumes a cosmology of 
D M = 0.3, fi A = 0.7 and H = 70 km s _1 Mpc~'. All errors are 
quoted at 68% confidence levels. 



2 CLUSTER SAMPLE & OBSERVATIONS 

In an attempt to decouple the effects of sample selection from pos- 
sible redshift evolution in X-ray properties, our targets were cho- 
sen to match our X-ray selected CNOC comparison sample in 
both redshift and mass. Velocity dispersion was used as a mass 
proxy for two-thirds of our targets, and initial masses for the re- 
maining third were estimated from cluster richness (Yee & Lopez- 
Cruz 1999; Gladders & Yee 2005). The ranges of redshifts and 
velocity dispersions spanned by our sample are 0.16 < z < 0.42 
and 709 < cr < 1390 km s -1 , compared to CNOC ranges of 
0.17 < z < 0.55 and 575 < a < 1330 km s _1 , Each target was 
observed in the X-ray by either Chandra or Suzaku. Table 1 lists 
the clusters in our sample along with their redshifts, the telescope 
used, obsids, and exposure times. 



2.1 Chandra Initial Processing 

Seven of our nine targets were observed with Chandra's ACIS-S 
(Advanced CCD Imaging Spectrometer) array in VFAINT mode 
during the period 2009 November 24 - 2010 August 3. The field 
of view of the ACIS-S aimpoint chip is 8'x8'. Chandra has excep- 
tional 0.5"spatial resolution, enabling very detailed investigations 
of cluster cores. The spectral resolution of ACIS is 120-130 eV over 
the energy range of interest to this study (0.3-7 keV). 

Preliminary reduction of the Chandra data was performed as 
in Hicks et al. (2008). After this initial cleaning, 0.3-7.0 keV and 
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Figure 1. Smoothed Flux Images. Adaptively smoothed 0.3-7.0 keV Chandra flux and 0.2-12.0 keV Suzaku counts images of our sample. RCS1 102-0319, 
RCS1 102-0340, and the second RCS2347-3535 panel are Suzaku observations. Circles denote calculated values of R2500 for each cluster. In each image, north 
is up and east is to the left. 
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Table 2. Cluster Positions and Detection Details 
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RCS2347-3535 
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23 47 49.2'' 
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11.6 
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a All positions are given for equinox J2000 

h 0.3-7.0 keV band, within R < 500 Kjl kpc 

c Suz.ak.it observation, R = 260" 

d Using Chandra data 

e Using Suzaku data 



0.3-2.5 keV images, instrument maps, and exposure maps were cre- 
ated for each dataset using the CIAO 4.2 tool MERGE_ALL and 
CALDB version 4.3.0. Point source detection was performed by 
running the tools WTRANSFORM and WRECON on the 0.3-7.0 
keV flux images. Data with energies below 0.3 keV and above 7.0 
keV were excluded due to uncertainties in the ACIS calibration and 
background contamination, respectively. These products were used 
solely for imaging analysis; our spectral analysis excludes all data 
below 0.6 keV. 



2.2 Suzaku Initial Processing 

The Suzaku X-ray Imaging Spectrometer (XIS) observed three of 
the clusters in our sample between May 31 and December 17, 
2008 including one that was later observed by Chandra (RCS2347- 
3535). Suzaku's field of view is 18'xl8', and it has a spatial resolu- 
tion of only ~ 2'. One of the advantages of Suzaku, however, is its 
excellent low-energy spectral resolution (60 eV at 0.25 keV). 

Suzaku data were reprocessed to include the most recent cal- 
ibration files using the HEASOFT 6.9 Suzaku software Version 
16 tool AEPIPELINE. Calibration sources were removed from the 
data using XSELECT, and data from XIS0, XIS1, and XIS3 were 
combined to produce integrated images in the 0.2-12.0 keV energy 
range. 



3 IMAGES AND SURFACE BRIGHTNESS PROFILE 
MODELLING 

3.1 Images and cluster centers 

Figure 1 contains smoothed 0.3-7.0 keV Chandra flux and 0.2- 
12.0 keV Suzaku counts images of each of the clusters in our sam- 
ple (produced by the CIAO tool CSMOOTH). Using these images 
we determined the location of the X-ray emission peak of each 
cluster. Optical positions were taken from RCS cluster catalogs, 
which were computed from the smoothed galaxy distribution used 
for cluster identification (see Gladders & Yee 2005, for details). 
Suzaku observations were included in this exercise despite having 
~ 2' spatial resolution. All clusters were found within 27" of their 
optical positions, and all Chandra observed clusters were found 



within 19", with a mean offset of 10.4". Table 2 lists optical po- 
sitions, X-ray positions, net counts (C - B) and signal-to-noise ra- 
tios derived from the method described in Appendix A. All clusters 
are detected at high significance. The median number of source 
counts per observation is 2911, with values ranging between 339 
(RCS 1615+3057) and 14,022 (RCS 1447+0828). 

3.2 Surface Brightness Profiles 

Radial surface brightness profiles were extracted in circular annuli 
from both targeted observations and blank sky background files us- 
ing 0.3-2.5 keV Chandra counts images. The cluster profiles were 
binned to ensure S/N > 3 in each bin and background-subtracted. 
None of our surface brightness bins are narrower than 1", therefore 
we do not expect to see any effects from Chandra's 0.5" FWHM 
PSF on our profiles. 

The surface brightness profiles were fitted to a beta-model: 

/(r) = / ^l + -j , (1) 

where I is the normalization and r c is the core radius. We compared 
the measured surface brightness (with error bars computed from 
photon statistics) to the average model surface brightness integrated 
within each radial bin while performing a Levenberg-Marquardt 
least-squares fit to the model, as implemented in the IDL routine 
mpcurvefit 1 . 

Best fitting model parameters for the seven C/janrfra-observed 
clusters are given in Table 3, and images of these fits are shown in 
Figure B 1 . Though many of the clusters exhibit some substructure, 
most were reasonably well fitted by the /? model (see Table 3 for 
goodness of fit data). /3 values for our sample are generally lower 
than the typical value of 2/3 (Neumann & Arnaud 1999), with a 
sample average of only 0.42 ± 0.01, suggesting that the gas is less 
centrally-concentrated in these targets compared to those studied 
by Neumann & Arnaud (1999). 

RCS2347-3535 contains a region of excess emission in its 
central 10-20 kpc, suggesting that it may harbour a small cool core. 

1 Available courtesy Craig B. Markwardt, 

http://cow.physics.wisc.edu/ craigm/idl/idl.html 



© 0000 RAS, MNRAS 000. 000-000 



X-ray Properties of Optically Selected Clusters 5 



Table 4. Spectral Fits (A = 2500) 
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a Suzaku observation, R = 260" 

^ Central 0.15 R500 removed from spectrum. 

c Best simultaneous fit to Chandra and Suzaku spectra. 

d Chandra normalization; Suzaku normalization for this fit is 1 1.5+2 * X 10 -4 cm 
e Using only Chandra data. 
f Using only Suzaku data. 



Table 3. /?-Model Fits 
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RCS0222+0144 
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None of the other clusters in our sample exhibit a significant central 
surface brightness excess above the /3 model. 



4 SPECTROSCOPY 

4.1 Initial spectral fits and redshift estimates 

Initial spectra were extracted from each point-source-removed 
event file in a circular region with a 300 h^J kpc radius (Chandra) 
or 260" radius (Suzaku). These values were chosen to provide an 
initial spectrum with as close to an R2500 extraction radius as possi- 
ble. In the case of Suzaku observations, the smallest recommended 
radius for extracting spectral data is 260". Individual spectra were 
extracted from each Suzaku XIS CCD and fitted simultaneously. 

All spectra were analyzed in XSPEC (Arnaud 1996), us- 
ing weighted response matrices (RMFs) and effective area files 
(ARFs) generated with the CIAO tool SPECEXTRACT (Chan- 
dra), and XISRMFGEN and XISSIMARFGEN (Suzaku). For the 
Chandra observations, blank sky backgrounds were renormalized 
to match the 9.5-12.0 keV count rates of each exposure, then ex- 
tracted in identical regions as source spectra. Suzaku backgrounds 
were compiled from the remaining data after calibration sources, 
point sources, and the target cluster were removed. 

Spectra were fitted with single temperature spectral models, 
inclusive of foreground absorption, using Cash statistics (Cash 
1979). Each spectrum was fitted with the absorbing column frozen 
at its measured value (Dickey & Lockman 1990). Metal abundances 
were allowed to vary. Data with energies below 0.6 keV and above 



7.0 keV were excluded from the Chandra fits, while Suzaku spectra 
were fitted in the 0.5-8.0 keV energy range (e.g., Bautz et al. 2007). 

Redshifts were fitted for the three clusters in our sample that 
do not have spectroscopic measurements. We were unable to con- 
strain a redshift for RCS 1330+3043, and therefore use its photo- 
metric redshift (z = 0.27) in further analysis. The fits of the other 
two clusters resulted in z = 0.334 ± 0.007 for RCS1 102-0319, and 
z = 0.160+°™ for RCS2150-0442, within 7% and 20% (respec- 
tively) of photometric redshift estimates obtained using the color 
of the red sequence (e.g., Gladders & Yee 2005). Fits with redshift 
fixed at these values were used in subsequent analysis (Table 4). 
The small uncertainties in these values 10%) do not substan- 
tially affect our analysis. 

A core-removed (0.15 R500) spectrum was extracted and fitted 
for RCS 1447+0828, indicating that it possesses a strong cool core, 
with an integrated temperature of 6. 8+5J" 3 , keV and core-removed 
temperature Tx=13+J keV (Table 4). Due to the central excess in 
its surface brightness profile, a core-excised spectrum was also fit- 
ted for RCS2347-3535; however, no conclusive evidence for soft 
central emission was found in this case. 



4.2 Scaled aperture estimate 

With an average of only 3000 counts per target, we do not have 
sufficient signal to undertake a full hydrostatic mass analysis, or 
even fit core-excised spectra out to Rsoo- Due to the shallow nature 
of our observations, we choose R2500 as the fiducial radius for our 
measurements. 

We employ scaling relations in order to estimate R^soo- The 
present sample contains objects in a variety of dynamical states 
and thus it is essential to use an appropriate scaling relation. We 
measure R2500 f° r our Chandra observations using Y x (= M g T x ) 
as a mass proxy. The decision to use Yx was motivated by two 
principal factors. Firstly, the simulations of Kravtsov et al. (2006) 
indicate that Y x is a robust mass proxy even in the presence of sig- 
nificant dynamical activity. Secondly, use of Y x facilitates compar- 
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Table 5. Sample quantities estimated within R2500- 
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fl Suzaku observation, Lx,2500 estimated within aperture of radius R = 260". 
° Core-excised temperature used to determine . 



ison of our optically selected sample to the X-ray selected REX- 
CESS sample, since REXCESS uses Y x exclusively. 

We used an iterative process to estimate the cluster masses. 
First we calculated central densities using the /? model parameters 
from Sect. 3.2 and initial spectral fits described in Sect. 4.1: 



4^ 2 , g (l+z) 2 



K 10' 



' c ,o 



0.82 Ami El 



(2) 



(Ettori, Tozzi & Rosati 2003). Here K is the normalization of the 
XSPEC model and El is the emission integral, estimated by inte- 
grating the (spherical) emission from the source out to 10 Mpc (see 
Ettori, Tozzi & Rosati 2003, for details). Values of n c0 were then 
used to determine total gas mass by integrating the cluster density 
profile out to R500: 



Pgas(r) = po 



1 + 



-3/J/2 



(3) 



where po = n c ,o fie nip, and \i e = 1.17. We then use the M-Y x 
relation of Arnaud, Pointecouteau, & Pratt (2007) to estimate the 
mass and radius of each object. The Arnaud, Pointecouteau, & Pratt 
(2007) M-Yx relation uses T[0.15-0.75 R 50 o] as well as M g [R 500 ]. 
We translate our T[R25oo] measurements to T[0. 15-0.75 R500] us- 
ing the mean T[R 250 o]/T[0. 15-0.75 R 500 ] ratio from the REX- 
CESS sample. R500 was then converted to R2500 using the sim- 
ple relationship R25oo= 0.44 R 50 o (Arnaud, Pointecouteau, & Pratt 
2005) 2 . The process of extracting spectra, fitting temperatures, de- 
termining gas mass, and re-estimating the mass and radius R2500 
was repeated until the radius values converged. 

Given the lower spatial resolution of Suzaku, in these cases we 
instead determined R2500 using the M-T x relationship of Arnaud, 
Pointecouteau, & Pratt (2005). The cluster RCS2347-3535 has 
been observed by both Chandra and Suzaku. Here Chandra spatial 
information was used in tandem with combined temperature fitting 
of both Chandra and Suzaku spectra. To confirm that the larger 
extraction region of Suzaku was not impacting best fitting tempera- 
tures, Chandra and Suzaku spectra were also fitted individually. All 
values of Tx (including the core-excised value) are consistent with 
one another within lcr errors. 



Use of a constant R2500/R500 ratio assumes that a constant concentration. 
This is a fair assumption given the mass range of interest (see e.g., Pointe- 
couteau et al. 2005). 



4.3 Scaled quantities 

Table 5 lists total gravitating masses (M lot ), gas masses (M g ), gas 
mass fractions (f gas ), and temperatures (T x ) within R25oo- Total 
gravitating masses estimated from the scaling relationships M- 
T x and M-Y x were found to be in statistical agreement with one 
another, with a slight tendency for the latter to be smaller. 

Bolometric (unabsorbed 0.01-100 keV) luminosities 
L x [R25oo] were calculated by extrapolating the best fit spectrum 
within XSPEC. Error bars on L x take into account uncertainties 
in R2500, temperature, abundance, and spectral normalization. 
To investigate whether different methods of determining R2500 
have a significant effect on resulting L x -T x relationships, we also 
calculate R2500 for the Chandra clusters using M-T x (Arnaud, 
Pointecouteau, & Pratt 2005), and re-extract spectra at those 
radii to obtain additional luminosities. We find that even when 
R2500 values differ by 20%, measured luminosities do not change 
appreciably (Table 5). Since we do not expect T x to change 
drastically on ~ 50 kpc scales at radii R J?R25oo , we find that 
L x -T x relationships are robust to different methods of determining 

R2500- 



5 CLUSTER SAMPLE COMPARISONS 
5.1 Context and comparison samples 

In our previous work (Hicks et al. 2008) we found that high- 
redshift, optically selected clusters had noticeably different X-ray 
properties than our moderate-redshift, X-ray selected comparison 
sample (CNOC). With our current moderate-z optically selected 
sample we are in a position to confirm whether these differences 
are due primarily to selection bias or to redshift evolution in cluster 
properties. We use the following comparison samples: 

• Ten high-redshift (0.6 < z < 1.0) optically selected RCS clus- 
ters selected from among the optically-richest of the 6,483 candi- 
dates detected in the first 90 deg 2 of the RCS (Gladders & Yee 
2005). The Chandra X-ray observations of these clusters are pre- 
sented in Hicks et al. (2008). 

• Fourteen moderate-redshift (0.2 < z < 0.55) X-ray selected 
CNOC clusters. This sample (Yee, Ellingson, & Carlberg 1996) 
is derived from detections in the Einstein Observatory Extended 
Medium-Sensitivity Survey (EMSS Gioia et al. 1990). CNOC 
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was originally chosen as a comparison sample because it has been 
extensively observed in optical, with galaxy redshifts of ~ 1200 
cluster members as well as detailed photometric catalogues avail- 
able (e.g., Yee et al. 1996). It is not, however, necessarily rep- 
resentative of X-ray selected samples as a whole. As the highest- 
luminosity clusters from the X-ray flux-limited, wide-area EMSS 
survey, the CNOC sample is almost certainly drawn from the ex- 
treme X-ray bright tail of the cluster distribution. The Chandra 
X-ray observations of these clusters are detailed in Hicks et al. 
(2006). 

• Thirty-one clusters from the R EXCESS sample (Bohringer 
et al. 2007). REXCESS provides homogenous coverage of the 
luminosity range 0.4 - 20 x 10 44 h~% erg s -1 in the 0.1-2.4 keV 
band (T x > 2 keV) over the redshift range 0.055 < z < 0.183. 
It is representative of the X-ray selected population, clusters hav- 
ing been selected in X-ray luminosity only, without regard for any 
other characteristic (apart from angular size, so that they would fit 
into the field of view of XMM-Newton). The X-ray scaling prop- 
erties of the REXCESS sample within R 50 o are discussed in Pratt 
et al. (2009). For the present paper, we have recalculated the X- 
ray temperatures and bolometric luminosities within R2500 defined 
as in Sect. 4.2 above. The resulting R2500, T x , L x , f gas > values are 
listed in Table 6. We also list the cool core (cc) and morphological 
classifications (Pratt et al. 2009). Objects are classied as cool core 
if central density E(zY 2 n e0 > 4 x 10~ 2 cm 4 . Objects are classi- 
fied as morphologically disturbed if the centroid shift parameter 3 
(w) > 0.01^500- The remaining systems are neither cool core nor 
morphologically disturbed. 

(Bohringer et al. 2010) 



5.2 Density Profiles 

Using the yS-model fits given in Table 3 and central densities found 
in Table 4, radial density profiles were produced, scaled by Rsoo, 
and compared to those of the REXCESS sample. The left hand 
panel of Figure 2 shows a comparison between the scaled density 
profiles of the moderate-redshift RCS and REXCESS samples. On 
average the scaled profiles of the RCS clusters are suppressed in the 
central regions compared to REXCESS. Note that this is a strong 
effect, reaching out well beyond R2soo- RCS 1447 is included in the 
calculation of the mean profile. Excluding it from the calculation 
would further suppress the average scaled density profile of the 
RCS systems in the central regions. 

The middle panel of Figure 2 shows the RCS clusters com- 
pared to the non-cool-core REXCESS subset only. Once again, 
there is a clear offset in the average density profile that extends out 
well beyond R95oo- 

The right hand panel of Fig. 4 shows the RCS clusters com- 
pared to the morphologically disturbed REXCESS subsample. 
Here the agreement between mean scaled density profiles is ac- 
tually quite remarkable. Note that there are two cool-core systems 
in the REXCESS morphologically disturbed subset, just as there 
is a cool core in the RCS moderate-z sample (RCS 1447), thus both 
samples span the whole range of cluster properties. Also, the mass 



3 The centroid shift parameter is denned as the standard deviation of the 
projected separations between the X-ray peak and the centroid at each ra- 
dius in the [0.1-1] R500 region (see e.g., Pratt et al. 2009; Bohringer et al. 
2010) 



ranges of the two samples match quite well, and this good agree- 
ment in mass is reflected in the similar temperature range. There- 
fore this effect is not likely to be the result of variations in gas mass 
fraction with mass. 



5.3 The L x - T x Relationship 

Studying the relationships between global cluster properties (L x , 
Tx, M tot , etc.) over a broad range in redshift allows us to inves- 
tigate the influence of non-gravitational processes on cluster for- 
mation and evolution. In addition, these relationships give us in- 
sight regarding cluster dynamical state and composition, as well 
as provide ready methods of comparison between different clus- 
ter samples. In this work we focus primarily on X-ray temperature 
(T x ) and luminosity (L x ), since they are tied most closely to the 
actual data and require fewer assumptions than extrapolated prop- 
erties such as total cluster mass. To facilitate comparisons between 
optically selected and X-ray selected clusters, we make use of our 
previous Chandra analyses of the CNOC (Hicks et al. 2006) and 
high-z RCS samples (Hicks et al. 2008). 

In all L x -T x relationships, L x has been scaled by the cosmo- 

r , -il/2 

logical factor E z = H(z)/H = K2,„(l + z) + Q.a\ and the rela- 
tionship has been fitted with the form 

log 10 (E:'L X ) = Cl+C2 log 10 T x , (4) 

with T x in units of 5 keV and L x in units of 10 44 erg s -1 . New 
best fitting relationships are determined using both the bisector and 
orthogonal modifications of the BCES algorithm in Akritas & Ber- 
shady (1996), with figures displaying bisector fits for ease of com- 
parison to our previous work. 

It is immediately clear from visual inspection of Figure 3a that 
our current targets are overall less luminous for a given temperature 
than the CNOC sample, in keeping with our previously-observed 
RCS sample. This result suggests that global discrepancies between 
the X-ray properties of X-ray and optically selected cluster samples 
are wholly attributable to selection bias. 

After correcting for self-similar evolution we find no convinc- 
ing evidence of redshift evolution in cluster properties within the 
complete X-ray observed RCS sample (0.15 < z < 1.0). To further 
explore this point, we fitted the RCS moderate-z and high-z samples 
together with the function 

log 10 (L x ) = A log 10 (EJ1.3) + B log 10 (T x /5) + C (5) 

The best fitting values are (bootstrap uncertainties): A = 0.43± 
1.15; B = 2.69 ± 0.84; and C = 3.23 ± 1.24, showing that we are 
quantitatively unable to constrain evolution with the current sam- 
ple. 

Figure 3b shows RCS and REXCESS L x -T x scaling relation- 
ships, for the entire REXCESS sample and for just the non-cool- 
core (ncc) clusters. For the REXCESS sample, cool cores define 
the upper envelope of the relation while morphologically disturbed 
systems define the lower envelope. The vast majority of the disper- 
sion comes from the cool cores. While REXCESS cool-core clus- 
ters look much like the high-L x CNOC sample, the non-cool-core 
clusters are situated similarly to the optically selected sample. 

It is clear that the normalization, in particular, that one will 
obtain for any L x -T x fit will depend strongly on the sample com- 
position. Note also that the intrinsic dispersion about the REX- 
CESS L x -T x relation is considerably larger than the 70% dis- 
persion about the L x -T x relation for quantities measured within 
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Table 6. R EXCESS X-ray Properties at R 2: 
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Figure 2. Density Profiles Moderate redshift RCS (dashed) vs REXCESS (dotted) density profiles, including mean density profiles for each sample 
(solid lines). Figure 2a: RCS vs full REXCESS sample. The density profiles of the RCS clusters are on average flatter in their central regions than the 
REXCESS systems, but are in good agreement with the X-ray systems at large radii (~R5oo)- Figure 2b: RCS vs non-cool-core REXCESS systems only, 
also exhibiting a slight offset. Figure 2c: RCS vs REXCESS morphologically disturbed subsample only. Here the mean profiles of the two samples are 
essentially identical. 



R500 (Pratt et al. 2009). The larger dispersion we see is due to mea- 
suring quantities in the core regions, where the variation in density 
is strongest. 

Building on the results of Section 5.2, we compare REX- 
CESS morphologically disturbed clusters to the RCS sample (Fig- 
ure 3c). Here we find an even closer similarity in X-ray properties, 
to the extent that even best-fitting relationships lie nearly on top 



of one another. This result is expected given the agreement in gas 
density profiles between our targets and the REXCESS morpho- 
logically disturbed subsample. L x -T x fits for a variety of sample 
combinations are reported in Table 7. 
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Figure 3. Lx-Tx Relationships. X-ray temperatures are plotted against cosmologically consistent, intrinsic bolometric luminosities within R2500- m a U panels 
diamonds represent high-z RCS clusters (z avg = 0.80), triangles denote the current sample (z avg = 0.30), and the dashed line shows the Lx-Tx relationship 
for the entire X-ray observed RCS sample. Figure 3a: Squares designate moderate redshift CNOC clusters (z aV g = 0.32), and the solid line shows their best 
fitting relationship. As in Hicks et al. (2008), the X-ray selected sample exhibits a higher normalization than optically selected clusters. Figure 3b: Circles 
represent the REXCESS sample. Here the solid line shows the fit to non-cool-core REXCESS clusters (filled circles) and the dotted line is the best fit to 
the entire REXCESS sample. Open circles denote cool-core clusters. Figure 3c: Here circles show only the disturbed REXCESS clusters, and the solid 
line represents their best fitting Lx-Tx relationship, which is notably similar to that of the RCS combined sample (dashed line). 



Table 7. Scaling Relationship: log 10 (E^Lx) = CI + C2 log 10 T x 

Sample Bisector Orthogonal 



Ci Co C\ C2 



CNOC" 


0.74 


+ 


0.08 


2.3 ± 


0.3 














RCS (high-z) a 


0.36 


+ 


0.06 


2.1 ± 


0.3 














RCS (total) 


0.34 
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0.06 


2.7 ± 


0.5 


0.34 


+ 
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2.i 


!± 
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0.1 


3.i 
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0.3 


0.47 


± 
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3.4 ± 


0.4 


REXCESS (total) 


0.63 


+ 


0.06 


3.1 ± 


0.3 


0.68 


± 


0.07 


3.', 


) + 


0.5 



a Hicks etal. (2008) 
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5.4 Gas Mass Fractions 

Using our estimated masses (Table 5) we calculate core gas mass 
fractions within R25oo- Our sample mean is .058 + 0.008 using Yx to 
determine total mass (error computed by dividing the standard de- 
viation by yN). In comparison, R EXCESS clusters have a mean 
gas mass fraction of 0.081 ± 0.003 for the complete sample of 31 
clusters, and 0.066 ± 0.004 when including only the 12 disturbed 
clusters; further evidence of similarities between optically selected 
clusters and morphologically disturbed X-ray selected clusters. 

We perform statistical tests comparing the distributions of 
f gas values for a subset of the clusters in our moderate redshift RCS 
sample with two subsets of the REXCESS sample. In order to 
compare clusters of similar masses, and to be consistent with Hicks 
et al. (2008), we restricted the temperature range in all 3 samples 
to be between 3.5-8.0 keV, resulting in 5 clusters from the RCS, 
17 clusters in the REXCESS sample, and 7 clusters in the REX- 
CESS/disturbed sample. The mean gas fraction and standard devi- 
ation from each of the subsamples are 0.050 ±0.005, 0.084 ±0.004, 
and 0.071 ± 0.005, for the RCS subset, REXCESS subset, and 
REXCESS/disturbed subset respectively. The RCS subset shows 
lower f sas than the REXCESS cluster subset, but similar f sas to the 
REXCESS/disturbed cluster subset, also visible in the histograms 
shown in Figure 4. 

We also performed a two-sided Kolmogorov-Smirnov test (the 
contributed KSTWO IDL routine, based on Press et al. 1992, 
Numerical Recipes), pairing the RCS subset with each of the 
REXCESS and REXCESS/disturbed subsets. The KS test com- 
paring RCS clusters with the complete REXCESS sample 
showed a higher probability of having been randomly drawn 
from different parent populations than the test involving only the 
REXCESS/disturbed subsample; however, such a comparison is 
severely limited by the small number of clusters. The KS probabil- 
ity comparing the 5 RCS clusters with the 17 REXCESS clusters 
was 0.013, while the KS probability comparing the RCS clusters 
with 7 REXCESS/disturbed clusters was 0.091. 



6 SUMMARY AND DISCUSSION 

We have presented a detailed X-ray investigation of a sample of 
moderate-redshift, optically selected clusters of galaxies from the 
Red-sequence Cluster Survey (Table 1). All of our targets were de- 
tected by Chandra and/or Suzaku at S/N>11 (Table 2) and were 
found within 27" of their optical positions. Chandra imaging re- 
veals that most of these clusters possess at least some degree of 
substructure (Figure 1). 

Surface brightness profiles were extracted for the seven 
C/zamira-observed clusters in S/N>3 annular bins. These profiles 
were reasonably well fitted by single /J-models (Table 3 and Fig- 
ure Bl). Cluster emission was modeled with XSPEC, beginning 
with a spectral extraction region of 300 kpc radius. The results 
of single-temperature spectral fits (Table 4), combined with gas 
masses obtained using best-fit /3-models, were used to determine 
Yx and consequently R2soo- This process was carried out iteratively 
until extraction regions and R2500 estimates were in agreement. 

Using both the M-Y x (Arnaud, Pointecouteau, & Pratt 2007) 
and M-T x relationship (Arnaud, Pointecouteau, & Pratt 2005), to- 
tal masses were estimated out to Risoo- While both relationships are 
overall consistent within our errors, we find M-Y x mass estimates 
to generally be lower than those obtained via M-T x for the clus- 
ters in our sample, not entirely surprising given that M-Yx uses gas 



mass information as well as temperature. Reassuringly, uncertain- 
ties in R2500 do not lead to large uncertainties in Lx[R25oo]- There- 
fore for X-ray observations for which only global T x and L x mea- 
surements are possible due to poor spatial resolution and/or low net 
counts, using M-Tx to determine R2500 results in reasonably accu- 
rate L x -T x relationships. Values of R2500, total mass, gas mass, gas 
mass fraction, and integrated luminosity are given in Table 5. 

The scaled density profiles (Figure 2) and global quantities 
of the moderate-redshift RCS sample are most consistent with the 
morphologically disturbed subset of the REXCESS X-ray selected 
sample. These density profiles explain trends in L x -T x relationships 
(e.g., Figure 3; Hicks et al. 2008) without having to resort to non- 
standard evolutionary effects. This result suggests that much of the 
non-standard evolution found by previous work may be due to the 
varying fractions of cool cores in different samples, which in turn 
are affected by the various methods of cluster selection. 

The properties of the scaled density profiles are reflected in the 
global quantities T x and especially L x , and vice versa. Low L x sys- 
tems have centrally-suppressed density profiles; high L x systems 
have centrally-peaked density profiles. Trends in the Lx-Tx rela- 
tion measured in an aperture corresponding to R2500 are similar to 
those measured at R 50 o (Pratt et al. 2009). cool-core systems pop- 
ulate the high-Lx side of the relation; morphologically disturbed 
systems populate the low-L x side. 

The fraction of cool-core clusters in the low-redshift RCS 
sample appears to be lower than in the REXCESS sample (1/7 
versus 10/31). This lower fraction is consistent with the differences 
seen in density profiles and scaling relations. However, the small 
numbers and potential mismatch in mass distributions keep this dis- 
crepancy in cool-core fraction from being statistically significant. 
If there is an increase in the cool-core fraction towards the present 
day, there is a possibility that the percentage of cool cores found in 
X-ray and optically selected samples might converge at higher red- 
shift. However, there are very few minimally-biased X-ray selected 
cluster samples that probe the redshift ranges we are interested in. 
The eROSITA survey (Predehl et al. 201 1) will definitely shed light 
on this issue. 

Overall we find that optically selected RCS cluster proper- 
ties span the entire range of those of massive clusters selected 
via other methods, but contain a higher fraction of objects with 
gas density profiles and global properties similar to those of the 
morphologically-disturbed systems in REXCESS. This result sug- 
gests that optical and X-ray selection do not sample exactly the 
same population of clusters. Recent results from SZ surveys such 
as Planck (Planck Collaboration et al. 2011), further suggest that 
X-ray selection may preferentially pick up centrally-concentrated 
systems. Selection effects such as these have the potential to affect 
the interpretation of cluster surveys that intend to use the evolution 
of the cluster population as a proxy for cosmic evolution. More in- 
vestigations are clearly necessary. 
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drawn from the same distribution. Right: Using only the 7 disturbed REXCESS clusters with 3.5 <Tx< 8 keV. 
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APPENDIX A: SIGNAL-TO-NOISE RATIOS 

To estimate the significance of RCS cluster detections in the X-ray, 
we used the same simple statistics as in Hicks et al. (2008). For 
the Chandra observed clusters, counts were extracted in the 0.3-7.0 
keV band from a 500 kpc radius region around the X-ray peak 
(C) and also from a region away from the aimpoint on the same 
chip which served as a background (B). Due to the ~ 2 arcminute 
spatial resolution of Suzaku, all analysis was performed using a 260 
arcsecond extraction radius, as recommended by the Suzaku Data 
Reduction Guide 4 . Therefore to determine S/N for our Suzaku ob- 
servations, counts were extracted in the 0.2-12.0 keV band from a 
260" radius region centered on the cluster (C), and the remainder of 
the chip minus calibration sources was used as the background (B). 
Obvious point sources were removed from each region in all obser- 
vations. Final signal-to-noise ratios were calculated based on divid- 
ing net counts, N = C - B, by the standard deviation, cr = VC + B. 
Using this method, all clusters were solidly detected at a signal-to- 
noise ratio greater than 1 1 (Table 2), with an average S/N of 41 for 
the sample. 



APPENDIX B: B MODEL FITS 



4 http://heasarc.nasa.gov/docs/suzaku/analysis/abc/ 
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Figure Bl. Surface Brightness Profiles. Background-subtracted radial surface brighness profiles for the 0.3-2.5 keV band accumulated in S/N> 3 annular 
bins for seven of the nine clusters in our sample. A solid line traces the best fitting single /? model of each cluster. Vertical dashed lines indicate R2500- Many 
of the profiles exhibit some substructure; however, most were reasonably well fitted by a standard fi model (see Table 3 for goodness of fit data). The surface 
brightness profile of RCS2347-3535 suggests that it harbors a modest cool core in the inner 10-20 kpc. 
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